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We theoretically study and experimentally report the propagation of nonlinear hysteretic torsional pulses in a 
vertical granular chain made of cm-scale, self-hanged magnetic beads. As predicted by contact mechanics, the 
torsional coupling between two beads is found nonlinear hysteretic. This results in a nonlinear pulse distortion 
essentially different from the distortion predicted by classical nonlinearities, and in a complex dynamic response 
depending on the history of the wave particle angular velocity. Both are consistent with the predictions of 
purely hysteretic nonlinear elasticity and the Preisach-Mayergoyz hysteresis model, providing the opportunity 
to study the phenomenon of nonlinear dynamic hysteresis in the absence of other type of material nonlinearities. 
The proposed configuration reveals a plethora of interesting phenomena including giant amplitude-dependent 
attenuation, short term memory as well as dispersive properties. Thus, it could find interesting applications in 
nonlinear wave control devices such as strong amplitude-dependent filters. 


Nonlinear dynamic hysteresis is involved in a wide 
range of acoustic effects recently observed in complex 
solids, often called mesoscopic solids (H. It is shown 
to originate from clapping and friction phenomena, in¬ 
duced by acoustic waves at the internal contacts and 
inter-grain boundaries of polycrystalline and Earth ma¬ 
terials EMU, damaged solids (30 and "model" granu¬ 
lar media El ED. This hysteretic nonlinearity enriches 
the phenomenology of the classical nonlinearity (small 
expansion terms of the smooth nonlinear stress-strain re¬ 
lation and geometric nonlinearity) with phenomena like 
nonlinear attenuation 120, memory CD and slow dy¬ 
namics Ifl2l[l3ll . 

One of the first observed manifestations of mechani¬ 
cal hysteresis is the shear coupling between two elastic 
spheres in contact lfl~4l IT5ll . Interestingly, periodic line- 
assemblies of spheres, also called granular chains, have 
attracted a growing interest in the last years for their 
rich nonlinear dynamics phenomena (harmonic gener¬ 
ation, solitons, breathers...) Ifl6l420l and the potentially 
wide applications in wave control devices ED- How¬ 
ever, up to now, dynamic hysteresis has not been studied 
in granular chains. In addition, all the media exhibit¬ 
ing dynamic hysteresis (rocks, disordered granular me¬ 
dia, concrete, composites...) have been also exhibiting 
at the same time other types of nonlinearities, classical 
or nonclassical. Consequently, although some model¬ 
ing has been proposed (22, [23j, it has not been possi¬ 
ble yet to observe pulse wave distortion by purely hys¬ 
teretic nonlinearity. Manifestations of hysteretic nonlin¬ 
earity were observed only in narrow frequency band ex¬ 
periments with sine waves, via for instance the shift of 
vibration resonances mm but never in the distortion 
of pulsed acoustic signals. The transformation of pulse 



Figure 1. (a) Experimental configuration for the torsional wave 
propagation (pure rotation of the beads around the axis z of the 
chain), (b) Signal acceleration at bead 15 versus frequency for 
a 15-bead chain, (c) Dispersion relation retrieved from the p 
resonance frequencies of the chain. 


profiles in media with classical quadratic nonlinearity, 
leading to weak shock front formation of the particle ve¬ 
locity profile, is one of the most classical observations 
in nonlinear acoustics EH 25 ]. Interestingly, the pulse 
distortion in media with cubic elastic nonlinearity, was 
observed for the first time only recently l26l . 

In this work, we report for the first time the transfor¬ 
mation of pulse profile in a medium with pure hysteretic 
quadratic nonlinearity, essentially different from the dis¬ 
tortion by classical nonlinearities. We start with the ob¬ 
servation and characterization of torsional (or pure rota¬ 
tional) wave propagation in a granular chain composed 
of identical magnetic spheres. Due to the pure torsional 
coupling at the contacts excited in a rotational motion 





along the axis z of the chain, pure nonlinear hysteretic 
behavior is observed. The latter is quantitatively char¬ 
acterized for a single contact in a resonance experiment, 
and then used for the modeling of torsional wave pulse 
distortion by quadratic hysteresis. A quantitative com¬ 
parison with the experimental distorted pulse profiles in 
a chain of 70 beads, shows the validity of the developed 
modeling. Extremely large nonlinear self-action (includ¬ 
ing self-attenuation and pulse deceleration) is demon¬ 
strated for parts of the pulse profile, depending on the 
previous loading by the pulse itself (a short term mem¬ 
ory). 

Experimental configuration. The medium is com¬ 
posed of 15 to 70 magnetic beads of diameter d = 2R = 
19 mm and mass m = 27 g (see FigJTJa)) 127. 28|. The 
magnetic poles of the beads are aligned along the chain 
(z-axis) and provide a magnetic force To = 54 N be¬ 
tween beads, normal to the contacts. F 0 exceeds largely 
the gravity force and thus the contact forces along the 
chain are considered identical in this study. The chain 
is excited at the top by a shaker (B&K 4810) coupled 
to the side of the bead n = 0 , which produces a rota¬ 
tional motion relative to the 2 -axis of the chain. This 
dynamic rotation is transmitted along the chain thanks to 
the torsional rigidity of the contacts, leading to torsional 
wave propagation. For detecting the waves, accelerome¬ 
ters (type PCB 352C23) are glued to the side of several 
beads and angular particle accelerations are recorded via 
an oscilloscope or a spectrum analyzer. 

The dynamics of the chain can be described through 
the following system of equations of rotational motion 
for each bead n interacting with its nearest neighbors n — 
1 and n + 1 through torsional moments, 


jd 2 Pn 
dt 2 


— M n — M n _ 1, 


( 1 ) 


where /3 n is the rotation angle of bead n, J = 
2mR 2 /5 ~ 9.85.10 -7 kg.m 2 is the moment of inertia 
of the beads, t the time, and M n the nonlinear torsional 
moment depending on the relative rotation angle between 
beads n+1 and n. Following the approach in 111411291 , the 
moment-angle relationship for oscillatory motion (a sim¬ 
ple periodic driving) between two spheres in contact can 
be approximated by a linear and a quadratic hysteretic 
functions (see the Supplementary Material), 


M n = M^ n + M‘ 
Kt S 'Ipn h 


( 2 ) 


Vni’n + ~ (ipl ~ Vn) 


where fi) n = fi n +i — fin is the relative rotation angle 
between two adjacent beads and is its magnitude, 


K t = d( 1 — u)Fo is the linear torsional constant, h = 
Eafi / (1 + z^/iTo is defined as the parameter of quadratic 
hysteretic nonlinearity, with \i the a priori unknown fric¬ 
tion coefficient of the bead’s material, E = 160 GPa 
the Young’s modulus, v ~ 0.24 the Poisson ratio of the 
bead’s material, and a = (3FqR( 1 — v^/AE) 1 / 3 ~ 
130 pm the contact radius. 

In Fig. [3b), the angular particle acceleration A 15 of 
bead n = 15, is presented versus frequency in the case of 
a 15 bead chain. In this finite-length chain, resonances 
are observed up to the cut-off frequency f c ~ 275 Hz. 
Simulating the linear dynamics of the chain with springs 
having a complex constant K t = 0.78(1 + 0.027i) N.m 
to fit the experimental losses, we obtain the theoretical 
curve in Fig.[ljb). The theoretical cut-off frequency f c = 
y/Kt/J / 7 r = 283 Hz is found in good agreement with 
the experimental one. From the experimental resonance 
frequencies f p (with p = 1,..., 15) of the finite chain, 
the dispersion relation can be retrieved. At the p-order 
resonance of the chain, the wavenumber is k p = (2 p — 
1)tt/AL, where L is the length of the chain. The relation 
f p — k p is evenly derived for the simulation results and 
the two dispersion curves are compared successfully in 
Fig.[ljc). These observations provide an estimate of the 
long-wavelength torsional wave velocity in such system 
c = dy/Kt/J — 17 m/s. 



Figure 2. Magnitude of the acceleration transfer function (de¬ 
tected acceleration A\ at bead 1 over the acceleration Ao at 
bead 0) of a "2 beads-1 contact" system for different excitation 
levels, experiment (a) and theory (b). (c) Relative resonance 
frequency shift as a function of the detected resonance ampli¬ 
tude (normalized to the one at the maximum excitation level), 
(d) Setup for the single contact characterization in torsion. 

In order to verify in our configuration the validity of 
the nonlinear moment-angle relationship ([ 2 ]), and to es¬ 
timate the parameter h, we characterize a single con¬ 
tact between two spheres in a nonlinear resonance ex¬ 
periment for rotational motion as depicted in Fig. |2jd). 
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The nonlinear resonance method has been widely im¬ 
plemented for characterizing the hysteretic behavior be¬ 
cause the classical quadratic nonlinearity does not con¬ 
tribute to the amplitude dependent resonance shift at the 
leading order I9lfl3l. The downshift of the resonance fre¬ 
quency with a linear dependence on the oscillation am¬ 
plitude can be attributed to hysteretic quadratic nonlin¬ 
earity mm, while the shift proportional to the square of 
the amplitude at lowest excitation levels - either to cubic 
elastic nonlinearity 130- 32] or to hysteretic nonlinearity 
described by Preisach-Arrhenius model l33ll . 

In Fig. [2ja), the experimental ratio between acceler¬ 
ation signals at bead 0 and bead 1 is shown as a func¬ 
tion of normalized frequency D = ///o, where fo = 
159 Hz is the experimental resonance frequency at the 
smallest excitation amplitude. Downward resonance fre¬ 
quency shift is observed as well as a nonlinear attenu¬ 
ation process EH. The relative resonance frequency 
shift A///o = ( f r — /o)//o is shown to scale linearly 
with the detected resonance amplitude in Fig.|2jc), which 
is consistent with quadratic hysteresis. For comparison, 
the theoretical transfer function obtained with the har¬ 
monic balance method for the moment-angle relation¬ 
ship Eq. ([2]), and neglecting contributions from higher 
harmonics 11341 !35i , is plotted in Fig.|2jb). A quantita¬ 
tive agreement between the theoretical and experimental 
resonance curves and resonance frequency shifts is ob¬ 
tained for/i^ 258±16 (which corresponds to a realistic 
friction coefficient /i ~ 0.32 =b 0.02). The fitted quality 
factor Qo = 37 is also in agreement with the imaginary 
part of K t used above. 



Figure 3. Experimental angular particle velocity 0 V = 0 v /0o 
at bead 1 for the smallest excitation amplitude and its approx¬ 
imation by a linear piece-wise function. T n are the positions 
of the signal zeros in time normalized by to = 5.77 ms so that 
Ti — 1, r n are the extrema positions in time, and S n are the 
extrema amplitudes normalized by 6o so that S\ = 1. 

After having extracted the linear properties of the 
medium and tested the quadratic hysteretic behavior for 


a single contact, we turn to the nonlinear propagation of 
torsional wave pulses. In order to avoid simultaneous de¬ 
tection of forward and backward propagating waves, and 
to mimic a semi-infinite chain, we use a 70-bead long 
chain. The detector is placed at the bead n — 18, and 
the wave propagates through 104 beads before coming 
back to the receiver after reflection by the free boundary. 
The pulse central frequency should be sufficiently low 
to avoid strong dispersion effects occurring at frequen¬ 
cies close to the Bragg frequency f c ~ 275 Hz, but also 
sufficiently high to limit its spatial extent and be able to 
distinguish incident from reflected pulses. Therefore a 
central frequency of 100 Hz has been chosen. The first 
bead velocity is shown in Fig. [3j corresponding to the 
time-integrated accelerometer signal. 

We describe the pulse distortion starting from the fol¬ 
lowing evolution equation, valid for any type of signal in 
a one-dimensional dispersionless medium where linear 
dissipation is neglected 122, 23, 35], 

de v i dM h de v _ 

2 d§ v dr { j 

Here, in accordance with the considered problem, 9 V = 
0 v /0q is the normalized angular particle velocity, 0q is 
the angular velocity amplitude of the first phase of the 
emitted wave packet at z = 0 , r = (t — z/c) /to where 
to = 5.77 ms is the duration of the first phase of the 
emitted signal, = Mh/(dK t ) and £ = z/z n i with 
z n i = c 2 to/hdOo the characteristic nonlinear length (37). 

The important term dMh/dQ v , having the physical 
sense of a normalized amplitude-dependent modulus 
(22ll23ll36) . and containing the hysteretic behavior, has 
been analytically derived from the Preisach-Mayergoyz 
(PM) model of hysteresis CD [38 39 ]. This phenomeno¬ 
logical model, developed initially in magnetism and then 
adapted to elastic waves, considers a large number of 
hysteretic elements (hysterons), each having two possi¬ 
ble stress states. The transitions from one state to another 
take place at two characteristic strains, one for each sign 
of the strain rate. An hysteretic medium can be repre¬ 
sented via a distribution of hysterons in the PM plane 
formed by two axes whose coordinates are the two char¬ 
acteristic strains of the hysterons. Unlike existing ana¬ 
lytical formulas for the quadratic hysteresis, this model 
has the ability to model the instantaneous memory stored 
in the hysterons’ states under the action of an arbitrary 
varying acoustic loading and thus, is applicable to arbi¬ 
trary signals such as the pulses observed here. 

The evolution equation ^ is then modified into a sys¬ 
tem of equations, derived for the linear piece-wise ap¬ 
proximated signal ofFig.0 The derivation using the PM 
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model of hysteresis is detailed in the Supplementary Ma¬ 
terials, and leads to the following system, 

dT n = S n d£/2, (4) 

dr n = (T n - r n )A n d£/(T n - T n ^ t ), (5) 
dS n = -A n S n d£/(T n - T„_i), (6) 

with A n = S n if S n S n —i and A n = ( S n 4- S n —^/2 

if S n ± Sri —i* The system 0-© describes respec¬ 
tively the small shifts of the positions in time of the ze¬ 
ros (dT n ), the extrema of the signal (dr n ), as well as the 
change in extremal values (dS n ) when a small change in 
dt; occurs, i.e. either a small distance or a small excita¬ 
tion amplitude change. So, by increasing step by step the 
amplitude £, the system |4])-([6]) describes the distortion 
by the hysteretic nonlinearity of the linear piece-wise ap¬ 
proximated pulse. We note here that the evolution of 
the parameters of the signal phase number n depends 
in general on phases n and n — 1 but not on previous 
(n — 2, n — 3...). Equations 0 and 0 have a clear 
physical sense. The local time delays are proportional 
to wave amplitudes as it could be expected for quadratic 
hysteretic nonlinearity (T|. Equation ([6]) describes local 
nonlinear absorption with an absorption coefficient pro¬ 
portional to wave amplitude and to the local "frequency" 
(inverse duration of the pulse phase), also theoretically 
expected for quadratic hysteretic nonlinearity Ifm In- 
teresting compression/stretching effects could take place 
for slowly modulated wave packets containing many car¬ 
rier periods. 

In Fig. |4ja) we plot the experimental angular veloc¬ 
ity signals for different excitation amplitudes at bead 1 
(z = 0 so by definition for £ = 0) and at bead 18 
(z = 1 8d = 0.323 m), i.e. for different values of £. 
For the smallest excitation amplitude, the characteristic 
angular velocity is #o — 0.17 rad/s and the experimen¬ 
tally found speed of sound is c = 20 m/s. As observed 
and predicted for waves in resonance experiments [ 2 - 
sm both effects of nonlinear softening and nonlin¬ 
ear attenuation are manifested here for the multi-phase 
pulses through stretching in time and normalized am¬ 
plitude decrease with increasing excitation level, respec¬ 
tively, Fig.[4|c,d). Additionally, it is striking to observe 
that different parts of the signal are not attenuated in the 
same way. For instance, the third extremum is attenu¬ 
ated faster than the second phase, being initially smaller 
in amplitude (see Fig.[l](c)), which is counter-intuitive. 
This effect is explained by the model in Eqs. ([4][6]) where 
the evolutions of signal phases amplitudes are coupled 
with the phase durations. For h = 800 ±160, all these 
observations are quantitatively captured by the modeled 
signal distortion from Eqs. 0-0. This value of param¬ 
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Figure 4. (a) Experimental angular velocity signals for differ¬ 
ent values of £ : at bead 1 (£ = 0 dashed line) and at bead 
18 for different excitation amplitudes (£ =0.5, 1, 1.3, 1.8 and 
2.4 for h = 800). (b) Corresponding theoretical angular veloc¬ 
ity signals, (c) Experimental (symbols) and theoretical (lines) 
amplitudes of the four first signal extrema versus £. (d) Exper¬ 
imental (symbols) and theoretical (lines) durations of the three 
first signal phases versus £. Error bars limits represent simu¬ 
lated results for h = 640 and h = 960. 


eter of hysteretic nonlinearity is different from the one 
evaluated above for the one contact - one bead system 
resonance h ~ 258, although of the same order of magni¬ 
tude. This deviation could be attributed to the differences 
in the setup (chain of 2 beads versus 70 beads) through 
the role of the weight, but also to the possibly different 
contact wears and their expected important influence on 
h through the friction coefficient (i. 

Conclusions. Existence and propagation of pure tor¬ 
sional waves in a granular magnetic chain have been re¬ 
ported. The torsional waves show dispersive properties 
associated to the periodicity of the medium. Through a 
single contact nonlinear resonance experiment, we veri¬ 
fied a nonlinear quadratic hysteretic relationship for the 
moment-angle (equivalent to stress-strain). Following 
this observation, we developed a set of equations based 
on the Preisach-Mayergoyz model of hysteresis to de- 
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scribe the torsional pulse propagation. The model’s 
results compare very well with experimental distorted 
pulses observed for increasing excitation amplitude in 
a 70 bead long chain. The reported nonlinear trans¬ 
formation of pulse profile in a medium with hysteretic 
quadratic nonlinearity, essentially extends the histori¬ 
cally observed distortion by quadratic nonlinearity m 
and the more recently observed distortion by cubic non¬ 
linearity [26]. 

In particular, we found a giant nonlinear attenuation 
(or shift in time) for each phase of the pulse, which de¬ 
pends on the characteristics (amplitude, duration...) of 
the particular phase of the signal and sometimes on the 
relative amplitudes of the previous phases. This signifies 
a short term wave-memory effect. The observed giant 
nonlinear attenuation effects and the short term memory 
could become key components of nonlinear elastic wave 
control devices and logic elements ESI28BD. 

The presented configuration could be used as a model 
medium to study fundamental nonlinear wave pro¬ 
cesses, e.g. frequency mixing, interaction of counter- 
propagating waves, self-modulation instability, still un¬ 
explored for hysteretic nonlinearity. At the same time, 
it could also be used to study novel nonlinear processes 
which do not exist with classical nonlinearities, like the 
nonlinear memory and the influence of the preparation 
of the initial hysteron distribution on the nonlinear pulse 
propagation. 
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